Contradicting effective mass scalings within the Skyrme energy density functional 

method 
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The problem of the effective mass scaling in the single particle spectra calculated within the 
Skyrme energy density functional (EDF) method is studied. It is demonstrated that for specific 
pairs of orbitals the commonly anticipated isoscalar effective mass (m*) scaling of the single-particle 
level splittings is almost canceled by an implicit m*-scaling due to other parameters in the Skyrme 
EDF. This holds in particular for an indirect m*-scaling of the two-body spin-orbit strength making 
the theory essentially unpredictable with respect to single particle energies. It is argued that this 
unphysical property of the Skyrme EDF is a mere consequence of the strategies and datasets used to 
fit these functionals. The inclusion of certain single-particle spin-orbit splittings to fit the two-body 
spin-orbit and the tensor interaction strengths reinstates the conventional m*-scaling and improves 
the performance of the Skyrme EDF. 
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Fundamental excitations in atomic nuclei are often 
characterized in terms of single particle and collective 
excitations. A successful nuclear theory is required to 
account for both kinds at a qualitative and quantita- 
tive level. The effective nuclear energy density function- 
als (EDF) irrespectively of their variants including local 
(Skyrme), non-local (Gogny) or relativistic mean- field 
(RMF) have developed toward high accuracy in recent 
years, both with respect to bulk properties and predic- 
tions of single particle properties [lj. Still, on a quanti- 
tative level, properties of single particle states are often 
better described by means of simple potentials like the 
Nilsson or Woods-Saxon. This is highly unsatisfactory, 
since one expects a full fledged EDF to describe nuclear 
properties at a level of accuracy superior to that of simple 
potentials. The uncertainty of present EDF with respect 
to nuclear properties has its roots in the fitting proce- 
dure of its parameters as well as the non local terms in 
the interaction. 

The EDF's are conventionally adjusted to reproduce 
ground-state bulk nuclear properties, see e.g. Ref.0. 
The datasets used to fit their parameters are dominated 
by nuclear data extrapolated to the thermodynamic limit 
and by nuclear binding energies in selected doubly magic 
nuclei. These fitting procedures are known to impair ba- 
sic buildings blocks of these theoretical models, in par- 
ticular single-particle (s.p.) energies. The physical rele- 
vance of s.p. energies provided by self-consistent mean- 
field (MF) approaches based on the EDFs or effective 
interactions is continuously contested. The debate has 
its roots in the non-locality of these approaches result- 
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ing in a low isoscalar effective mass — ~ 0.8 [3( which 
in turn scales the s.p. level density g in the vicinity 
of the Fermi energy ef according to the simple rule: 
g(sjr) — > ^t<7(£_f). This simple rule applies strictly to 
homogeneous nuclear matter only. In finite nuclei the ef- 
fective mass depends on r and the giep) — * ^prgi^F) scal- 
ing should be considered as an idealization. The effective 
mass scaling leads to a dichotomous and in fact highly 
uncomfortable situation. Indeed, in spite of the fact that 
it makes all applications of the self-consistent MF meth- 
ods to low-lying nuclear excitations rather dubious such 
calculations are carried out routinely and published often 
without even a word of comment or justification. 

The aim of this Letter is to demonstrate that the prob- 
lem of the effective mass scaling within the effective the- 
ory is far more intricate than anticipated. It turns out 
that the effective mass dependence of the calculated sin- 
gle particle spectrum can depend on the structure of the 
EDF, on the strategies employed when adjusting the in- 
teractions and upon the choice of dataset used in these 
fits. We will demonstrate first that conventional func- 
tionals which are, as discussed above, fitted almost ul- 
timately to bulk nuclear properties have an built in im- 
plicit effective mass scaling of certain coupling constants 
including in particular the two-body spin-orbit strength. 
For many modern parameterizations of the Skyrme force 
this mechanism is strong enough to cancel the effect of 
the s.p. level density scaling caused by the low effective 
mass with respect to the calculated specific particle-hole 
excitation. For example the — /V/2 splitting in A~44 
mass region [4| calculated using forces having effective 
masses ranging from 0.7 < — < I yield similar result in 
spite of the expected effective mass dependence. In the 
next step we will show that this counterintuitive result is 
a mere consequence of the fitting strategies and that by 
shifting the attention from bulk to single-particle prop- 
erties in the process of fitting of the nuclear EDF param- 
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eters one can both remove the artificial m*-scaling from 
the two-body spin-orbit strength and reinstate the con- 
ventional and anticipated m*-scaling in the calculated ph 
excitation energies. 

All the calculations performed in this Letter are based 
on Skyrme-force inspired local EDF (SEDF) of the form: 

«to= E (nr n (r)+H° dd (r)) , (1) 
t=o,i 

where 

nr D = C t p p 2 + Cf p p t Ap t + (2) 
Clp t T t + c/j? + C t VJ Pt V • J t , 

Hf d = C|s t 2 + C t As s t • As* + (3) 
C t T s t .T t + C t J j2 + C t Vj s t -(Vxj t ). 

It depends on isoscalar t = and isovector t = 1 time- 
even p t , Tt, and J t , and time-odd s t , T t , and jt, lo- 
cal densities for which we follow the convention intro- 
duced in Ref. see also Refs. d, 0] and references 
cited therein. All numerical results presented below were 
obtained using the HFODD code of Ref. [8j and the 
following set of the Skyrme forces: SLy4 L SLy5i [j, 
SIII L [§], SkM* L H3, SkXc L [10], MSkl L 0, SkP L 
and SkO_L [13]. The subscript L indicates here that the 
original time-odd functional coupling constants C s ,C As , 
and C T were replaced in the calculations by the coupling 
constants reproducing the empirical Landau parameters 
in accordance to the prescription given in Refs. pi flU- 

In spite of the fact that we will concentrate here on the 
implicit m*-scaling of certain functional coupling con- 
stants it should be mentioned that some of the SEDF 
parameters do scale or depend upon m* explicitly. The 
explicit effective mass scaling is well established for the 
(i) C s and C T functional parameters through the fit to 
the empirical Landau parameters [HHH and (ii) for the 
isovectorial coupling constants Cf and C{ through the 
fit to the empirical symmetry energy strength [lfij . 

An almost ideal playground to investigate the impact 
of fitting procedures on the performance of the effective 
forces with respect to the ph related observables is offered 
by fully stretched [/ 7 " 2 ]/ mol and [^3/2/7/2 k m „ states in 
Ny^Z nuclei in the A~44 mass region where n denotes 
number of valence particles in the fa n sub-shell pt] . In a 
scries of publications 0, [H M, El 

we have demonstrated 
that these states represent one of the best examples of 
almost unperturbed single-particle motion and that the 
energy difference: 

AE = E([f? /2 ] Imax ) - E([d-^+ 1 ] Imax ) (4) 

constitutes a very reliable probe of various properties of 
the EDF. In particular, it can be used to readjust the 
time-odd components C s , C As and C T to comply with 
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FIG. 1: The difference between theoretical and empirical val- 
ues of the energy differences AE defined in Eq. Q. Cal- 
culations have been done using different popular parameter- 
izations (see legend) of the Skyrme functionals with spin- 
dependent coupling constant readjusted to reproduce the em- 
pirical Landau parameters in accordance to Refs. [4l.ll5j. 



the empirical Landau parameters leading to an unifica- 
tion of the theoretical predictions for AE t h for various 
popular Skyrme paramctcrizations [4]. 

The calculated values of AE t h, see Fig.HJ show a strik- 
ing and completely unexpected feature. The mean devi- 
ation of the theoretical predictions with respect to the 
empirical values is very similar for such forces like SLy4, 
SLy5 ~ 0.7), SIII ~ 0.8), SkO ~ 0.9) and 
SkXc (^- ~ 1.0). In spite of the fact the the effective 
masses differ by as much as 30%, we obtain for these pa- 
rameterizations SE = AE t h — AE exp ~ — 550 keV with 
a rms deviation of a w 70 keV. This result is indeed ex- 
tremely puzzling since the anticipated influence of a naive 
77i*-scaling with respect to a ph excitation energy of or- 
der of AE exp w 5.5 MeV is, for the analyzed set of forces, 
estimated to be more than one order of magnitude larger, 
exceeding er w 1 MeV. 

Unexpectedly, the observable appears to be very 
robust. Indeed, the structural purity of the terminat- 
ing states in the A^44 mass region reveal a firmly es- 
tablished hierarchy of different physical contributions to 
that quantity and in turn allows to identify the physical 
source of the cancellation of the m*-scaling in the func- 
tional. In Ref. 3 we established the hierarchy of three 
different components, the energy scale hu, the strength 
of the spin-orbit (SO) and the £ 2 (surface) term using the 
Nilsson Hamiltonian. This schematic model clearly indi- 
cates the dominant influence of the mean SO potential 
on AE. This observation is neatly correlated with self- 
consistent SHF models, see the Appendix below where 
we give a numerical proof for the case of the Skyrme 
functionals. Indeed, the isoscalar strength of the Skyrme 
one-body SO potential emerging from the two-body SO 
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FIG. 2: The isoscalar strength of the two-body spin-orbit 
interaction plotted as a function of the isoscalar effective mass 
parameter for a few popular Skyrme force parameterizations 
studied here including SLy4, SLy5 SkO, SIII, SkXc, SkP, and 
SkM*. Stars mark the Skyrme forces having non-zero tensor 
terms. 



interaction Wo = —2C^ J appears to be almost perfectly 
linearly correlated with the effective mass m* parameter 
as shown in Fig. [2] Large values of Wq which character- 
ize low— to* forces tend to reduce the ^3/2-/7/2 splitting 
thus compensating the TO*-scaling effect on the N=Z=2Q 
magic gap. As shown in Ref. [201 ] a similar cancellation 
takes place for the 231/2-^5/2 and / '5/2-39/2 splittings in the 
A=16 and A— 80 mass regions, respectively. For these 
splittings the schematic Nilsson model indicates an in- 
disputable dominance of the SO term over the surface 
(or, more precisely, the ~ £ 2 ) contribution. On the other 
hand, for the P1/2-59/2 splitting in the A=80 mass region, 
the SO and the 1 2 contributions are predicted to be sim- 
ilar. For this particular case also the self-consistent SHF 
calculations provide results which are qualitatively differ- 
ent than in the previous cases indicating a more complex 
dependence. 

The implicit mass scaling of the two-body Skyrme 
spin-orbit strength has unexpected and serious conse- 
quences: Within a single theoretical framework two con- 
flicting scalings are present. Indeed, particle-hole exci- 
tations associated with the spin-orbit partners such as 
/V/2 — /5/2 scale directly with Wq. At the same time 
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ph excitations scale according to the 



Wq = —Wq as shown in Fig. [3] At first glance this puz- 
zling situation seem to have no satisfactory and unique 
solution within the conventional Skyrme EDF. Although 
a reduction in Wq (see open symbols in Fig. [3]) clearly 
improves the agreement to the data it can hardly be ac- 
cepted as a reasonable solution. Indeed, the figure in- 
dicates that the empirical Ae(fi/2 — /s/2) splitting (ex- 
perimental data are marked by thick horizontal lines in 
Fig. [3]) is reached first by reducing Wq strength in pa- 
rameterizations characterized by large-TO* values. In con- 
trast, agreement to the empirical Ae(/V/2 — (^3/2) splitting 
is obtained after reducing the Wo strength in 1ow-to* pa- 
rameterizations. Therefore, a conventional theory can 
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FIG. 3: Upper part shows the calculated energy splitting be- 
tween vfr/2 — vd-i/i single-particle states versus Wq calcu- 
lated using different Skyrme forces. The lower part shows the 
u f7/2 ~ u fs/2 spin-orbit splitting versus Wo calculated using 
the same set of the Skyrme forces. Black (open) dots mark 
results obtained using standard (10% reduced) strength of 
the two-body spin-orbit interaction, respectively. The figure 
clearly demonstrates the presence of two conflicting scalings 
within the Skyrme model. See text for further details. 



not reproduce simultaneously both empirical values. To 
state it differently: the standard set of Skyrme forces ap- 
pear to be incomplete when confronted to experimental 
data. 

Fig. [3] suggests that the simultaneous agreement for 
both Ae(/ 7/2 - d 3/2 ) and Ae(/ 7/2 - / 5/2 ) splittings 
calls for parameterizations having large to* and drasti- 
cally reduced Wo. Drastic reduction of Wo will spoil, 
however, the relatively good agreement for our high- 
spin observable AE, see Fig. [TJ Hence a compensa- 
tion mechanism is required. This mechanism exist, it 
is not new [2l|, [52] and it is provided by a tensor com- 
ponent. However, until very recently, there has been 
little need to invoke the strong tensor interaction in 
EDF's, see 0, [H, HE H, E3, M, HI. In our re- 
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FIG. 4: The energy difference of theoretical and empirical 
AE t h — AE exp energy splittings given by Eq. Q. Upper part 
shows the results obtained using SLy4^T, SkP^T, SIII.lt, and 
SkOLT parameterizations with modified spin fields and the 
two-body SO and tensor strengths readjusted to the empir- 
ical — f/5/2 splittings, see table [I] Lower port shows 
the differences AE* h — AE exp between the isoscalar effective 
mass scaled theoretical splitting AE* h = ^AE t h and the 
experimental value. See text for further details. 



cent contribution [301 ] we have shown that the tensor in- 
teraction can be rather unambiguously fitted using the 
iV/2 — fb/2 spin-orbit splittings in three key nuclei in- 
cluding: the isoscalar spin-saturated 40 Ca nucleus, the 
isoscalar spin-unsaturated 56 Ni nucleus, and the isovec- 
tor spin- unsaturated 48 Ca nucleus. Unlike in the other 
studies cited above our work firmly revealed the need 
for a simultaneous drastic reduction of the two-body SO 
strength. Such a procedure definitely changes the phi- 
losophy behind conventional fitting strategies by shifting 
the attention from mass dominated gross features to pro- 
cedures including carefully selected single-particle states 
which are used to adjust specific coupling constants in 
the functional. 

We are now in a position to verify the consistency 
of our new fitting strategy using high spin terminat- 
ing states. We have therefore repeated the calculations 
for the energy differences AE of Eq. fH but using pa- 
rameterizations SIIIlt, SkOLT, SLy4iT, SkP^T, and 
SkXc^T- Subscript LT indicates that these parameteri- 
zations have: (i) time-odd spin-fields readjusted to em- 
pirical Landau parameters and (ii) modified tensor and 
two-body SO strengths. New values of the tensor and 
the two-body SO functional coupling constants are col- 
lected in Tab. Q] All other coupling constants of these 
parameterizations are kept to their original values. 

The results of our calculations are depicted in the up- 
per part of Fig. |4] It is interesting to observe that all 
large-m* parameterizations including SkOLT ( — ~ 0.9) 
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TABLE I: Spin-orbit C V J and tensor isoscalar Co and isovec- 
tor C{ functional coupling constants adopted for different pa- 
rameterizations. These modified parameterizations are sub- 
sequently used in the calculations presented in Fig. [4] 
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in C A p (surface), C T (effective mass), and C VJ (spin-orbit) 
coupling constants in the Skyrme functional. 



and SkP^T ~ 1-0) are giving satisfactory agree- 

ment to the data. In the particular case of the SkXc^T 
(— ~ 1.0) force (not shown in the figure) the mean 
deviation from the data equals to 33keV only! At the 
same time low-m* parameterizations including SLy4^T 
(2£ ~ 0.7) and SIII LT ~ 0.8) strongly overestimate 
the data. This result clearly follows the conventional wis- 
dom related to the m*-scaling of the s. p. level density. To 
visualize this even better we have rescaled the theoretical 
energy difference by the effective mass AE^ h = — AEth 
and have plotted the difference AE^ h — AE exp in the 
lower panel of Fig. |4j Three out of four curves depicted in 
this figure follow closely each other and the experiment. 
The curve representing the SLy4iT is slightly below the 
trend, indicating too small a reduction of the two-body 
SO strength (by 35%) done in Ref. 30]. The quantity 
AE can therefore be used for further fine tuning of the 
SO strength and indeed, a reduction of the SO strength 
by 40% shifts the SLy4 curve by 250keV up as expected. 
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One has to remember however, that the entire effective 
mass scaling concept in finite nuclei is only an idealiza- 
tion. Last but not least let us observe that all curves 
shown in Fig. 0] do not reveal any clear isotopic and/or 
isotonic dependence, in contrast to the previous calcula- 
tions presented in Fig. [1] This indicates a clear improve- 
ment of the isovector channel due to the presence of a 
strong isovector tensor component in the functional. 

In summary, we have investigated the impact of the fit- 
ting procedures and the datasets used to fit the Skyrme 
energy density functionals with respect to their spectro- 
scopic properties. We have demonstrated that the use of 
parameterizations fitted to reproduce bulk nuclear data 
in the thermodynamic limit and to binding energies of 
selected double-magic finite nuclei to analyze spectro- 
scopic data may lead to rather poor results. This is 
due to an implicit internal m*-scaling of, in particular, 
the two-body SO strength which is a mere consequence 
of the fitting procedure. We have further demonstrated 
that conventional Skyrme forces having the two-body SO 
strength and tensor coupling constants fitted directly to 
the empirical s.p. — vjsii SO splittings behave ac- 

cording to the expected m*-scaling law for particle-hole 
excitations. The present study seems to confirm the com- 
mon expectation that effective interactions with large ef- 
fective mass have a superior performance for calculations 
of spectroscopic data. It sends however a clear message 
that this conclusion is strongly dependent on the fitting 
process for effective forces. Involving a larger set of single 
particle data, in particular selected high spin terminating 
states reveal the need for a considerable reduction of the 
conventional two-body spin-orbit term and at the same 



time the requirement for a strong tensor term introduced 
in Ref. 0. 
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I. APPENDIX 

In order to demonstrate the hierarchy of various con- 
tributions to the observable ([!]) in A^44 mass region we 
have performed additional calculations using the Skyrme 
EDF ©-©• The results are depicted in Fig. The 
upper part of the figure shows relative changes 6(AE) in 
the energy differences AE caused by ±5% changes in C Ap 
(surface), C T (= -&) (effective mass), and C VJ (= C v ^) 
(spin-orbit) coupling constants with respect to their orig- 
inal values. The lower part illustrates relative changes 
in the total binding energy SE([fy 2 ]i max ) of the aligned 
state [/™/ 2 ]/ maa . caused by these variations in the coupling 
constants. This figure clearly shows that the hierarchy 
established in Ref. 0| on the basis of the schematic Nils- 
son model is correct. The leading contribution to S(AE) 
is indeed the spin-orbit term. The contribution coming 
from the effective mass term is circa 2.5 times smaller and 
is clearly non-perturbative as it impacts binding energies 
by ^lOMeV. In the same spirit we have also investigated 
the influence of ±1% variations in the bulk energy pa- 
rameter C p on both AE and E([f™ /2 ] Ima J. These tiny 
variations in C p impact binding energies by ~15MeV 
and, at the same time, affect AE only by ~±30keV. 
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